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In a scenario with two almost mass degenerate heavy sterile Majorana neutrinos with a mass ~ 1 
GeV, we present the semileptonic decay widths of heavy charged pseudoscalars such as B mesons, 
either lepton-number-violating (B^ —>■ or lepton-number-conserving {B^ —>■ ^^e^7r^), 

mediated by such on-shell neutrinos. It had been explained in the literature that such decays may be 
eventually detected, and that they can present even detectable CP violation effects. In this work we 
point out that, in addition, such decays may present detectable effects of heavy neutrino oscillation, 
allowing us to extract the oscillation length and thus the heavy neutrino mass difference AMn, as 
well as a CP-violating Majorana phase. 

PACS numbers: 14.60St, ll.SOEr, 13.20Cz 


I. INTRODUCTION 

The neutrinos can be either Majorana or Dirac particles, although most of the neutrino scenarios suggest that the 
neutrinos are Majorana particles. In this work we will assume that they are Majorana. As a consequence, they can 
induce not just lepton number conserving (LNC) but also lepton number violating (LNV) processes. Such processes 
are neutrinoless double beta decays in nuclei [T], specific scattering processes HU and rare meson decays [THU]. 

For neutrinos with masses, neutrino oscillations were predicted some time ago [TB]. Oscillations of active (light) 
neutrinos were later observed [nHH], with the conclusion that the first three neutrinos have nonzero but light masses 
< 1 eV. The oscillations are sensitive only to mass differences, while neutrinoless double beta decays and rare meson 
decays can help with the determination of the absolute mass of the light Majorana neutrinos. The best present upper 
bounds on the absolute masses of the light neutrinos are obtained from cosmology rrii, >0.23 eV ^D] . 

The light neutrino masses can be produced via the seesaw mechanism |2Ij where more than three neutrino flavors are 
required and where all of them are Majorana. The light neutrinos in these seesaw scenarios have masses ~ AI|)/AI_r 
(< 1 eV), with A4jj being an electroweak scale or lower; the heavy neutrinos are very heavy, with masses Ain >> 1 
TeV, their mixing with active neutrino flavors being very suppressed ^ AAd/AAh (^ 1). Other seesaw scenarios exist 
where the heavy neutrinos have lower masses Mn < 1 TeV, Refs. [32], and even < 1 GeV [31 I33H37] : their mixing 
with the standard model flavors may be less suppressed than in the original scenarios. 

CP violation in the neutrino sector in scenarios with nearly degenerate heavy neutrino masses had been investigated 
in scattering processes in the literature |28j (resonant CP violation), as an effect coming from the interference of 
tree-level with one-loop effects from the neutrino propagators. Further, CP violation in the leptonic [T31 [15] and 
semileptonic meson decays [HHii] was investigated with a simpler, effectively tree-level, formalism and again in 
scenarios with nearly degenerate on-shell heavy neutrino masses, and where the decay width matrix of the massive 
neutrinos was assumed to be diagonal. 

Scenarios with nearly degenerate heavy neutrino masses (and CP violation effects) appear in various models, in 
particular in the neutrino minimal standard model (izMSM) |23ll29j where these neutrinos generate baryon asymmetry 
of the Universe while an additional lighter neutrino (with mass ~ 10^ keV) is responsible for the dark matter. Some 
general frameworks of low-scale seesaw (Ml El] with more than two heavy neutrinos can also explain the baryon 
asymmetry (but not the dark matter) and can have larger values of the heavy-light mixing than in z/MSM. In such 
models the almost mass degeneracy of Majorana neutrinos is preferred in the sense that it allows larger heavy-light 
mixings. 

In this work we discuss neutrino oscillations in semileptonic decays of heavy pseudoscalar mesons (such as B, Be, 
Dg) mediated by two on-shell Majorana neutrinos Nj {j = 1,2) which are almost mass degenerate. Similar effects 
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have been investigated recently in leptonic decays of such mesons, in Ref. [S], where a quantum field theoretical 
generalization of the Wigner-Weisskopf approach [32] was implemented and used. Our approach is simpler, and the 
results obtained are hopefully easier to interpret. 

In Sec. [IT] we present the results for the decay widths T and the effective decay widths Tea for the considered meson 


LNV semileptonic decay processes 




TT^ and their LNC counterparts 






the scenario with two on-shell almost mass degenerate heavy neutrinos Nj , 


all without the neutrino oscillation effects. 
Further, we also include the differential effective decay width dTeff{L)/dL for such processes, where L is the distance 
between the production {v-Nj) vertex and the decay vertex {Nj-e-ir). In Sec. Ill we then extend these expressions by 


including the oscillation effects of the on-shell neutrinos. In Appendix [A| we show the consistency of the oscillation 


amplitude method applied in Sec. |III| with the more usual quantum mechanics approach to oscillations. In Sec. |IV| we 
then estimate numerically the oscillation length and describe the conditions under which the oscillation modulation 
of the differential effective decay width dTee{L)/dL can be measured. In Sec. we indicate how the magnitudes 
\Bii\rj I of the heavy-light mixing parameters may be determined by such measurements. In Sec. 
results. 


VI 


we summarize our 


II. THE DECAY WIDTH EXPRESSION 


In this section we present formulas for the LNV and LNC semileptonic decays of charged B mesons, of the type 
B —>■ /ieTT, mediated by heavy sterile on-shell neutrinos. The formulas for the decay width of the LNV decays of this 
type, B^ —>• in the case one heavy neutrino, were presented in Ref. [10]. They were extended to the case of 

two almost degenerate heavy on-shell neutrinos in Ref. m (see also Ref. El): in the context of CP violation. For a 
review we refer to m- We will use these formulas, and also the formulas for the decay width of the LNC decays of 
this type, B^ —)• /i=*=e^7r=*=. 

There exist various scenarios with sterile neutrinos. Of particular interest is the z/MSM of Shaposhnikov et al., 
Refs. [531 [55]. This model contains two almost degenerate Majorana neutrinos Nj (j = 1,2) of mass ~ 1 GeV and 
another lighter neutrino lyx of mass ~ 10^ keV, as well as the three light neutrinos vj of mass < 1 eV. The striking 
advantage of this model is that it can explain simultaneously the existence of neutrino oscillations, dark matter 
and baryon asymmetry of the Universe. We also wish to point out that there exist more general (less constrained) 
frameworks of the low-scale seesaw, which explain baryon asymmetry but not the dark matter. In such models, larger 
values of the heavy-light mixing m are allowed than in i^MSM, and the case of almost mass degeneracy of Majorana 
neutrinos is preferred |31| since it allows larger mixings. 

The two flavor neutrinos and {i = e, /i, r) can be represented as 

3 

zzg = Bf.jVj + Bgpj^Ni + Bf,pf^N 2 + ■ ■ ■, (la) 

3 

t'/i = B^jVj + A^i + Bfj_]s[^N2 + .... (lb) 

1=1 


The coupling of these two flavor neutrinos to the corresponding charged leptons and has a part which contains 
the coupling to the heavy almost mass-degenerate neutrinos Ni and N 2 


^eWN — 2^^ L + BeN2^2) + h..c. — iFi + h.C. 

dlfiWN = (S^ATi A^i + Y2) -I-h.C. = +h.c., 

where we denoted by Afi and A /2 the e- and /r-flavor analogs of the heavy neutrino mass eigenfields Nj (j 


A/i = BiiNi + B12N2 

M2 = B21N1 B22N2 


Bak — 




-Rn 


A„ = v/|RaiP + |R«2P (a,A: = l,2), 


(2a) 

(2b) 


1 , 2 ) 


(3a) 

(3b) 


where in Bak the coefficients 0 = 1,2 stand for e, respectively; and fc = 1, 2 for Ai, AI 2 , respectively. The considered 
mechanisms for the LNV decays R^ —1 with on-shell Nj’s {j = 1,2), are those in Fig.[^ for the LNC decays, 

—1 /i^e^TT^ are those in Fig. Although the LNV decays have also the crossed channel (i.e., the ones where the 
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FIG. 1: The LNV decay B+ /i+e+ TT via exchange of an on-shell neutrino Nj [j = 1,2). 



FIG. 2: The LNG decay —>■ tt"*" via exchange of an on-shell neutrino Nj (j — 1, 2). 


vertices of ^ and e are exchanged), we assume here that the measurements can distinguish these two channels (since 
/i e), by reconstructing invariant masses from the detected final state particles. 

For these processes, we will consider the scenarios where the two heavy neutrinos and N 2 are almost mass- 
degenerate (AMat <C Mn = M]S!^) and are on shell. We will consider only the neutrino couplings ([^, with no 
components of the other mass eigenfields (thus no light mass eigenfields vi, V 2 , r'a), because we will assume that Ni 
and N 2 are the only neutrinos which are on shell in these processes. This is then reflected in our definition of “heavy” 
flavor states Nj, Eq. ( [3a] ). We stress that neutrinos which are off shell in these processes give, in relative terms, 
completely negligible contributions and will thus be ignored. 

In the more general case of the LNV decay —>■ ^ {£j = e, fi, r; M and M pseudoscalars), with 

i\ 7 ^ £2 and neutrino N on shell, the corresponding partial widths can be written as 


r(M± ^ £f A ^ £±£±m't) = E (4 ^ I 2 ), 


where 


F Mn ^i/ 2 (^ \ \l /2 I 1 y 2/^2 1 r^f '\ 

r = 04^2 ^ (i. 2 /Af> 2 /b) I 1 - —> — I Q[yN-,y£r,ye 2 ,y ) 


vn yN 


and the notations used in Eq. (j^ are 


Hyi , y 2 , y 3 ) = y? + yi + yf - 2 j/i2/2 - 21/22/3 - 22/32/1 


VN = 


Ml 

Ml 


M? 


yts = 


Mm 


y = 


MM 

Ml 


and the function (5(2/Ar;2/ti)2/f2i2/ ) is (cf- Ref. [T5] ) 


, (£s — £ 1 , £ 2 ) 


Q{yN-,yei,yi2,y') = f(2//v - yfi)(yiv - yG)(i - y/v - (1 - — +— 

yN yN 


+ [ - yGyG(i + y' + 2yiv - y^i -yM- yllyN - y') + yiM - yM 
+y^i(i + yiv)(y/v - y') - y^s]! - yiv)(yAr + y')] | • 


= 2 [(i“y^)y^+yA(i + 2 yiv-yG)] 


/ o y^2 (I \ 

VN-y - 2y/2-(y - y^J 

y/v 


(4) 

(5) 


( 6 a) 

( 6 b) 

( 6 c) 


(7a) 


(7b) 
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In Eq. ( 6 a), /m and Im' are the decay constants, and Vq^q^ and Vq^q^ are the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix elements of pseudoscalars and M ^. 

We notice that in the considered specific case (£i = and £2 = e) we have yt,^ « yi,^ « 0, and expression Q 
simplifies according to 


Q{yN-,yei,ye 2 ,y') ^ QiyN-,0,0,y') = ^[yN{i-yN)][yN -y']- 


The results 0-0 can be written in an equivalent form 


1 


ih ^ h), 


r(M± ^ efN ifefM =f) = —r(M± ^ £f fv)r(iv ^ £fM 

Tat 

where the widths of the two decays are 

TiM^^efN) = ^£fN), 

T{N^£fM'^) = \Bi,n\^T{N ^£^m'^), 

and the expressions for the corresponding canonical widths T (i.e., widths without the mixing factors) are 

T{M^ -)■ £fN) = X^/'^{l,yN,ye,) [(1 - yN)yN +yeA^ + 2yN - yij], 

1 -.2 nr 1 2 1 \l /2 ( 1 y y ^2 


( 8 ) 

(9) 

( 10 a) 

( 10 b) 


( 11 a) 


r(iV ^ £fM T) = —Glfl,,\Vq^qf— aV 2 1 , ^ [(m^ + m 2J(M^ - M^, + Ml) - 

iOTT ivi]^ \ y^v Vn / 




IOtt 


Vn yN 


yN-y' - 2ye2 - —ijj' - ye^) 

yN 


( 11 b) 


where again the notations 0 were used. We notice that the algebraic facto rizati on of the Q function, Eq. ( [7b| , yields 
the factorization ([^, as can be seen by inspection of the expressions ( 11 aI and ( 11 b). 

It can be checked that the result for the LNC processes —>■ £fN M ^ is the same as the result (0-( 10) 


r{M^ 


£fN^£f£fM'^) = —T{M^^£fN)T{N^£^M'^) 
Tat 


f'i'i 


( 12 a) 


\Be,N\ \Bi^n\ y{M+ ^ e+N)T{N ^ £+m'-) (£1 £ 2 ), (12b) 

Tat 


where the canonical decay widths (E’s) are again those of Eq. (11). 

We recall that we will consider the scenario with two on-shell neutrinos Nj (j = 1,2), and with almost degenerate 
masses: IAMat] Mat^, where AMj^ = Mjq^ — In this case, it turns out that the expression for the LNV decay 
width becomes more complicated, cf. Ref. [13]. With the notation ([^ for the mixing coefficients, it can be written in 
the following form: 


r(R=‘ 


= T{B'^ —>■ fi^N)T(N —)■ e“'' 7 r ) x { 
4 


\B^Ni \'^\BeNi I 

TaTi 


I Bfj,N 2 1 I Bgj^^ 

^N2 


Mr -lB^NA\BeNA\B^ 2 N 2 \\BeN 2 \ (s{y) cos ^ ^sind^^NV)^ ^^ 3 ^ 

(1 7 Vi + 1 Af 2 ) V y J 

where Mpf = Mpf^ « Mn 2 , the angle 02 i(I'Nv) is a combination of the phases of the heavy-light mixing coefficients 

gi(LNV) ^ arg(R^Ar2) + arg(ReJV2) - arg(R^Ari) - arg(ReAri)- (14) 

The functions S{y) and rj{y)/y appearing in Eq. in are functions of the parameter y = AMjq/T at only, where Tat is 
the arithmetic average of the total decay widths of Ni and ^ 2 ^ 


1 ^ ^ AMat 

= 7; i^Ni+BN2), y = —’ 

/ i at 


(15) 


^ For simplicity we assume that the (2 x 2) decay width matrix {TN)ij of the (near mass degenerate) neutrinos A^i and N 2 has no (or 
negligible) off-diagonal element [(rjv)i 2 ]- This assumption was also taken in Refs. |12H15| . 
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y y 


FIG. 3: The suppression factors ri{y)/y and S{y) as a function of y = AMat/F n'- (a) for 0 < 1 / < 1 (on the linear y scale); (b) for 
1 < y < 10 (on the logarithmic y scale). Figures taken from Refs. |13[ll5j . For j/ < 0 we have S(y) = S{—y) and ri{y) = v{—y)- 


and the functions are presented in Fig. [^for y > 0. It is straightforward to verify that the invariance of expression 
under the exchange of the roles of Ni and N 2 means that for y < 0 we have S{y) = S{—y) and r]{y) = r]{—y). 
The factors (functions) S{y) and 'rj{y)/y represent the effects of the iVi-fV 2 overlap in the decay width, in the real and 
imaginary parts of the N 1 -N 2 interference terms, respectively. Therefore, 5{y) and go to zero when |y| 3> 1, 

i.e., when no overlap. While both the 5{y) and the 'ri{y)/y function were obtained in Ref. |13j numerically, it can 
be argued that ri{y)/y is a simple function [15], r]{y)/y = y/(y^ + 1), and this agrees with the numerical results of 
Ref. [13]. For details, we refer to Refs. [T31 ITS] . 


Finally, the expressions for the total decay widths Fatj and Fappearing in Eq. (13) are 


Tat- = ICjTi\[{Mi\[) , 


(16) 


where we denote Mn = ~ Mm 2 and 




(17) 


and the factors Kj 


P U = 1)2) contain all the dependence on the heavy-light mixing factors 


fCj — Af zN 


BeNj P + \Bf2Nj P + A/'t AT \BrNj 


(18) 

neutrino 
Ref. [T3] 
neutrino 


In this expression, the coefficient functions = J\feN {(- = e, li, r) depend only on the mass Mat of the 

Nj-, these coefficient functions are ~ 10°-10^. The curves for AfmiMpf) as a function of Mn were presented in 
for the case of Majorana neutrinos, and in Ref. m for both cases of Dirac and Majorana neutrinos, in the 
mass interval 0.1 GeV < Mat < 6.3 GeV. 

As mentioned earlier, in addition to the above LNV decay width, there exists also the LNC decay width F(i3^ —>■ 
/i^e^TT^), which in the case of scenario of one on-shell neutrino N coincides with the LNV expression ([^. In the 
scenario with two on-shell almost degenerate neutrinos Nj , the expression is slightly different from the LNV equation 
(13); namely, only the angle [Eq. (14)] is now replaced by the following angle: 


^(LNC) ^ _ arg(ReAr2) - arg(R^Ari) + arg(ReAri). 


(19) 


We will consider, from now on, the case when there is an almost degeneracy of the two heavy neutrinos (jAMArj <C 
Mpf = and at the same time the degeneracy jAMArj is significantly larger than the (extremely small) decay 

width Fat: 


IAMatI Mat and |y| = 


IAMatI 

Tat 


( 20 ) 


In this case, we can see from Fig. that the functions 5{y) and 'r]{y)/y become very small. 


Therefore, the N 1 -N 2 
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overlap term in r(i3^ —>■ /i^e^Tr^) becomes negligible. As a result, Eq. (13) reduces to the following form: 

2 

r{B^ ^ ^ lt=^e=^7r=F) 

i=i 

= T{B+ ^ ^l+N)T{N ^ e+TT-) 


TaTi 


rATa 


(21a) 

(21b) 


The decay width presented hitherto does not contain an important suppression (acceptance) factor. Namely, the 
on-shell neutrino Nj travels before decaying. The decay will be detected if the on-shell neutrino decays during the 
passage of the neutrino through the detector. If the length of the detector is L, then the probability Pn of decay of 
N there is 


Pjv(T) = 1 - exp (- - — 

V tn^n 

'^/{tnJnPn) if Pv‘C 1 . 




(22a) 

(22b) 


In the second identity of Eq. ( |22a| ) we took into account that L = /S^t where /3n (^ 1) is the velocity of the N 
neutrino in the lab frame. Furthermore, 7Ar = (1 — is the Lorentz lab time dilation factor, typically jn > 2. 

The N lifetime in the rest is tn = I/Eat. Therefore, 'YnT'n is the N lifetime in the lab frame. The formula (22b) 
holds if Pn ^ 1, and we will assume this to be the case in the considered cases. 

This decay-within-the-detector probability Pn has been discussed in Refs. n HTHU [m I3M35] It is convenient to 
define the corresponding canonical, independent of mixing, probability Pn 


Pn{L) = 

, T{Mn) 

Im X - 

(23a) 

In 

^ Pn(.L) k. 

f L\ 

(23b) 


in Fig. as a 

function of Mn- 




FIG. 4: The canonical probability Pn, as defined in Eq. (231, as a function of the neutrino mass Mn, with the Lorentz lab 
time dilation factor chosen to be Jn [= (1 — ~ 2. 


The effective (true) decay widths and branching ratios are those multiplied by Pn- However, since we have two 
different (but almost mass degenerate) neutrinos Nj, we have for each of them a different decay probability 


PnAL) « ( X PnK:j = ^Tn^, 


Im J 


In 


(24) 
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where ICj is given in Eq. (18), and in the second equality we used the relations (23) and (16). The canonical 


probability Pn, Eq- (23), is common to both neutrinos Nj because they have practically the same mass and thus the 


same kinematics (and hence the same Lorentz factor jn)- The coefficients 10°-10^) in ICj are common to both 

neutrinos Nj (because they have a practically equal mass); but the mixings can be, in principle, quite different 

for the two neutrinos, and thus the two mixing factors Kj (j = 1, 2) may differ significantly from each other. 

Combining the probabilities (24) with the decay width ( |21b ) leads to the effective (true) decay width, where the 
dependence on the two decay widths Tjv cancels out: 


reff(B± ^ ^±e=^7r=F;L) « ^ ^ (25a) 

1=1 

« -^T{B+ ^ ^i+N)f{N ^ e+TT-)[\B^N,\^\BeN,? + \B^N,?\BeN,?]- (25b) 

INPN 


This implies that the effective differential decay, with respect to the distance L between the two vertices of the process, 
is 


—reff(B=^ ^ 


^ T{B+^ ji+N)V{N ^ e+n-)[\B^ 


InPn 




\Bii.N2?\BeN2?] 


(26) 


which is independent of the distance L. 

For the LNC processes B^ —i the result is the same as in the above LNV processes, due to the equal¬ 

ity of the LNC decay width ( [I^ with the LNV decay width in ([^ and (10) [cf. also Eq. ([II|)]. Therefore, when 
5{y), \rj{y)/y\ ^ 1 [i.e., when the conditions (20) hold], we have 


dL 


reff(B^ 


y^e^TT^\L) = —reff(B=‘= 
dL 

,±p=F^±. T,\ — 


reff(B=^ -1 y^e^TT^\L) = reff(B=^ 


Ai=^e=^7r=F;L), 


(27a) 

(27b) 


The present upper bounds for the p mixing coefficients appearing in these expressions, in the considered mass 
range Mjv « 1-5 GeV, are \BiMj\^ ^ 10“'^-10“"‘, cf. [8].^ 


III. THE EFFECTS OF NEUTRINO OSCILLATION 

In the previous section, important effects of neutrino oscillation of the propagating on-shell neutrino were not 
accounted for. As we will see, these effects lead to a modulation, i.e., the B-dependence of the effective decay widths 
obtained in the previous section, where L is the distance traveled by the on-shell neutrino between its production and 
detection points {L « P^t). 

We will follow the lines of the approach of Ref. [35] to neutrino oscillations. For the LNV decays B+ —)■ —)■ 

/i+e+TT” of Fig. the relevant interactions at the first (production) vertex are —B*^^JZ°7''(1 + ^^)NjW^'^\ and the 
neutrino state produced at this vertex is 

l^)(B+) B;^JiVi(m)) +B;^JiV2(p^J), (28) 

where the momenta of the two physical on-shell neutrinos are slightly different from each other, because jAMjvj yf 0 
(we recall that jAMiv] <C M^r). We have 

+ {p%^Y, (29) 

where the restriction to one spatial dimension (z) was made, because the processes with oscillation require the neutrino 
to propagate far from the production vertex. At the second vertex of the LNV process Fig. the relevant coupling 
is B*pf _Nj-j^{l — 75)eWj^^ The detection of the neutrino there can be described by an operator at the detector 
space-time location z = (t,0,0,B) where L « ^Nt- This operator is the annihilation operator B*j^_b(_is!.)(pMj'-iZ) = 


2 


The present upper bounds for are higher than that, but they are expected to become significantly lower in the future. 
















S:xp{—ipN^ ■ z) acting at the aforementioned component ^ ^(Arj)(piVj)^|0) (j = 1,2). Since 

b{N){PN)b{N){PNy \0} = constjO), this implies the following detection amplitude^: 


AiB+ p+e+TT ; L) ~ exp(-ip7Vi • z) + exp(-ipAr2 • z). 


(30) 


The L dependence of the effective (true) decay width of the considered process is proportional to the absolute square 
of the above amplitude 


d' j.,(osc) 

dZ 




-+-+--; L) = 


/i ' e ' TT 


rr\T\B+ 


^+e+7r"; L < L' < L + dL) ^ \-A{B+ 


+-+--)|2 (31a) 


II' e TT 


1 p(osc) 

dL 

f\BeNj 1 ^ + 2Re [B'^ff^B*pj^B^N2BeN2 exp [i{pN2 - Pm) ■ z]] | (31b) 

The superscript (osc) indicates that this is the (differential) effective decay width with oscillation effects included. 
The oscillation term, in comparison with expression (|26| , is new and introduces L-dependence in the otherwise L- 
independent differential decay width dTes/dL of Eq. ( |26[ ). This oscillation term comes from the interference term in 
the square of amplitude (301. Therefore, by comparing the obtained expression (31) with (26), we can obtain the 
complete expression for the effective differential decay width with oscillation effects included 

1 




r(B+ ^ ^+ 7 V)r( 7 V ^ e+TT") 


X \B^,N, ?\BeN, ? + 2Re [S*5eJV2 exp [i{pN 2 - Pn,) ■ | -(32) 

The oscillation term here contains two on-shell 4-momenta = (Ejv^, 0, 0,p^.) (j = 1,2) which are related by the 
on-shellness conditions pNj • PNj = and by the condition 


13 n2 — Pni = 


PN2 


^ 3\2 


Pm 


En2 


0 . 


(33) 


This condition comes from the following interpretation. The and N 2 amplitudes interfere at L if both of them 
are appreciable there. The neutrinos Ni and N 2 , in general, separate as they travel from their production to their 
detection vertex. Interference is then possible there only if this separation IAL12I = \{I3n 2 — /3Ni)\t (with; t « L/Pm) 
is smaller than the spread of the wave packet ALwp = cf. Ref. [36] 


|/3jv 2 - /3jVil ^ , ..X 

\1}N2+Pm\ t ^ 


(34) 


Stated otherwise, the following hierarchy is assumed: 




FIG. 5: Graphical representation of the hierarchy equation (351 of the lengths AL 12 , ALwp and the detector length L. The 
two interaction vertices (the production and the decay vertex of the neutrino N) are denoted as vi and V 2 , respectively. Note 
that at the production vertex (ui) the wave packets of Ni and N 2 are not mutually displaced, unlike in the decay vertex (U2). 


® We use the metric (1, —1, —1, —1) for the scalar products. It is our understanding that the authors of Ref. m use the metric (—1,1,1,1). 
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|ALi2 | ( = 


|^jV2 - Pni\L 

Pn 


< ALwp (=/3 jvAT) < L, 


(35) 


cf. also Fig. 

In order to express the oscillation phase 4)(L) = {pn 2 ~PNi) ■ z in Eq. (32) (34)] in a convenient form, condition 

(33) can be used. Since = P%^/Em 2 and = p%^/Eni are close in value, hence they are close to the value of 
Pn = {p %2 Ep%^)/{En 2 + EMl)- It can be checked to see that the latter velocity is practically equal to the arithmetic 
average {1/2){(3n2 + Pni)- Therefore, 


z = (t, 0,0, L) « t(l, 0,0, Pn) = 


{En 2 + Eni ) 


{pn2 +PNi) 


(36) 


Therefore, the oscillation phase is 


<(>(T) = {pn2 -PNi)- z = t- 


(M2 ) 


{En-i + Em2 ) 


tM 


N 


N- 


AM 

Em 


-^Mm^=L^, 
Pm Em PnIn 


(37) 


where it was taken into account that = M^_, and M^^ — = 2Mm^Mm (where AMm = Mm 2 — Mm^ and 

|AMv| «C Mmi = Mm)- As stressed in Ref. m. this expression for the oscillation angle is valid always, not just for 
relativistic neutrinos Nj, whenever relation (34) is fulfilled. For example, if the neutrinos are nonrelativistic, we have 
4>{L) « {L/Pm)AMm- The obtained oscillation phase allows us to define the oscillation length Lose as 

2ttPm1n 


0(Tosc) — 27r 


Using expression (37), the differential decay width 




InPn 


^ " AMm ■ 

can now be written in a more explicit form 

T{B+ p+N)T{N e+TT") 


(38) 




^\BsMf 


■ 2|R^Ari ll-SeiVi ll^eATa I COS ( L 


AMm , ^(LNV)^ 


(39) 


where the constant phase is defined in Eq. (|l4[). We can integrate the differential decay width (39) over the 


dL length to the full length L between the vertices. If L ^ |Tosc|) this then gives the full effective decay width of 
Eq. (25b) because the oscillation term ~ cos(^(L) + 6*12) gives a relatively negligible contribution when integrated 
over several “oscillation wavelengths” |Losc|- If, on the other hand, we do not assume L ^ |Tosc|, the integration of 
expression (39) gives 


F 


(osc) , „ + 


eff 


—>■ ; L) 


InPn 

Lose. 


ttL 


r(R+ ^ fi+N)r{N 

I BfiMi 11 BeNi 11 Bfj,M 2 11 BeN: 


e'^TT 




sin 


2,^ + e(“''>)-™(»S'‘''’) 

-^osc / 


• (40) 


Until now we considered the case of oscillation effects in LNV decays —J, /i+e+7r“. It can be checked that for the 

charge-conjugate LNV decays B~ —> p~e~'!T~^ the previous derivation can be repeated, with the only replacements 
Bg^, !-)• BiMj and B^Mj Bl^_. Instead of Eq. (30) we now have 


A{B p e 7r+; L) ~ B^^MiBeNi exp(-zpAr^ • z) + exp(-ipAr2 ’ z)- 


(41) 


This implies that in result (39) we now get 02i '—>■ —^*21, so that we can extend the results (39) and (40) to both LNV 
cases (i?±) 


_c(_j^(osc) / „± ± ± =F. T\ ^ ^ 


r(R+ ^ p+N)T{N e+TT") 




\BeNi\ +2|B^ArJ|5eAril|B^Ar2l|-BeAr2|cOS 27r-- ± 0. 


q(LNV)\ 


^21 




(42) 
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InPn 

Lrisr. 


ttL 


r(B+ ^ ^i+N)T{N 

I 11 ^eNi 11 11 BeN: 


e'^TT 


^ 1=1 

(LNV)\ 


sin 27r 


±e. 




(43) 


For the LNC processes B^ —i /i^e^7r=*= (cf. Fig. [^, in the case of no oscillation effects the results for the decay 
widths are the same as for the LNV processes; cf. Eqs. ([^-(27). When oscillations are accounted for, the results are 


almost the same as in the just considered LNV processes, except that for the decay amplitudes [cf. Eqs. (301 and (41) 
for LNV case] we have some of the heavy-light mixing elements B^m complex-conjugated and others not 


A{B+ /i+e" 7 r+;L) 
A{B~ —>■ fj,~e'^TT 


;L) 


B*^NiBeNi exp(-ipjv^ • z) + B*^^^BeN 2 exp(-ipjv2 ’ z), 
Bf^NiB*^^ exp(-zpAr^ • z) + exp(-ipjv2 • z). 


(44a) 

(44b) 


This then leads to the following results, in analogy with the LNC results (42 )-(|43|) where now only the phase angle 
021 ^^^ gets replaced by a different phase angle 02given in Eq. (19): 




—r(B+ 


InPn 


^i+N)T{N e+TT") 


^ \B^ 2 N,?\BeN,? + 2|i?^wJ|i?eA,J|i?^w,||i3e^,| cos ^ ± X45) 


p(osc)(^d 


/i^e^TT^; L) 


InPn 

LctfiC 


l^lBeNf 


ttL 


r{B+ fi+N)T{N e+TT") Iei^ 

2-7^±9ir')Tsin(«ir’)l|- (46) 

-^osc / J j 


fiNi I \-DeNi \\-DfiN 2 I \-^eN2 


All the formulas with oscillation effects, derived in this section, can be extended in a straightforward way to the 
oscillation effects in the semihadronic decays with two equal flavors of produced charged leptons, i.e., —)■ ^ 

and —>■ more specifically, B^ —>■ and B^ —>■ (cf. Sec. [v|. 

In Appendix we show that the wave function approach of Ref. m (cf. also [5S]) to the considered LNV and LNC 
processes with on-shell neutrinos is consistent, within their approximations, with the amplitude approach presented 
here and based on the method of Ref. [36] . 

A question may appear why the usual (^-matrix) approach, leading to the results of the Sec. [llj did not give 
the modulation (oscillation) effects obtained in this section. The x-coordinates of fields are integrated over in the 
^-matrix approach of Sec. |nj reflected by the use of initial and final states with definite momenta {Sp = 0). The 
uncertainty relation implies then 6x = 00 . Therefore, the location of the vertices remained undefined in the approach 
of the previous section. On the other hand, if the location of the vertices is to be determined in an experiment by 
precision Sx{= 6L) ^ 0.1 mm or better, then the corresponding precision in the determination of the momenta is 
Sp'> l/(5x ~ 10“^ eV. 


IV. OSCILLATION LENGTH AND MEASUREMENT OF THE MODULATION 


For the described oscillation modulation to be well defined and detectable, several conditions have to be fulfilled. 
Among them is the hierarchy (351 between the length L (the distance between the production and the decay vertices), 
the width L^p of the wave packet, and the separation AL12 between the two wave packets at the second vertex 
(cf. Fig. and Ref. [36]). Yet another necessary condition for the detection of the oscillation is that the maximal 
detected length L between the two vertices (we will call it simply the total detector length, Lmax = Adet) is larger 
than or comparable with the oscillation length |Aosc| [Eq. (38)]. For the measurement of the oscillation modulation 
effects in practice, the case [Aoscj ~ Adet is more convenient than jAoscj < Adet, i.e., 


/ 2t: \ 

\ [AMjvj / 


(— Adet). 


(47) 
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TABLE I: Presently known upper bounds for the sqnares \Bin\'^ of the heavy-light mixing matrix elements, for various specific 
values of Mjv- We excluded the upper bounds for B^nI^ from the 0vPj3 decay, which are uncertain due to possible cancellation 
effects. See also Figs. 3-5 of Ref. [39]. For each upper bound, the corresponding experiment (reference) is indicated. 


MArlGeU] 


IFivP 

Ib.atP 

1.0 

3 X 10-" (gOj) 

1 X 10-4 (im) 

3 X 10-4 (HD) 

2.1 

4 X 10-5 

3 X 10-5 (|43|) 

2 X 10-4 

3.0 

2 X 10-5f|l5)r 

2 X 10-V|T5)T 

4 X 10-5 (|42l) 

4.0-5.0 

1 X 10-4 IpSfY 

1X10-4 rpfr 

1X 10-4 r|i5)~ 


Furthermore, if the decay probability PNj{Ldet) for the decay of Nj (j = 1,2) within the detector [Eqs. ( [^ and 
is significant, i.e., if Pn {Ldet) ^ 1; then the oscillation is not well defined because it disappears within one 


©] i 

oFdess 


oFdess oscillation cycle due to the decay of Nj. Therefore, for the oscillation to be well defined, we have to require 


PNjiLdet) "C 1. This means, according to Eq. (22b) and using Eq. (47), the following: 


/I r I —N, 
U osci -) 


Pdet ^ 


PnIn 


Ni 


(48) 


This implies that we have l/|AMjv| «C I/Fat^. (j = 1,2), meaning that the condition |y|(= \AMn\/T]s[) » 1 of 
Eq. (20) is fulfilled when we have well-defined and detectable oscillation.^ We recall that this condition (|y| ^ 1) 


was assumed throughout the derivation of the oscillation formulas of the previous section so that the (otherwise 
problematic) overlap terms with S{y) and r]{y)/y factors in expression (13) could be neglected. 

The oscillation length can be estimated in the following way. Let us assume that the near mass degeneracy {y ^ 1) 
is in the interval: 1 ^ |?/|(= IAMtvI/Fat) < 10 ^, i.e.,^ 


IAMatI < IO^Ftv. 


(49) 


Furthermore, let us take that in the total decay widths Fat^ , Eqs. (16)-(18), the dominating contribution in the mixing 
factors JCj is from f'-component, i.e., ICj « (j = 1,2; £ = e or /i or r). Stated otherwise, we assume that 

is the largest among the mixings \Bf^Nj\‘^ and iBrNjl"^- Then, we have 


FAf. = 


J^lN 

IF 


f\B^ 

10-5 


f ■^in\ 

^ IF" 


IB, 


INi 


10-5 


X 4.57 X 10"^® GeV 

1 


43.5 m 


(50) 


For Mat = 1-5 GeV, and taking N to be Majorana neutrino, we have N^n ~ ~ 6-10, and A/)- at ~ 3-5 

(cf. Refs. [131H]); hence, the factor A/^jv/10 in Eq. (50) is ^ 1. The factor iR^Ar, p/10-5 in Eq. (50) can be ^ 1 
larger or smaller; cf. Table |T] for some present upper bounds. The oscillation length (38) can then W estimated 


or 


I Lose I - 


27r|gAf| ^ |r/v| 27r 

Mat IAMatI Mat lO^F, 
\pn\ 1 


10 


Mat 10 \^feN 


^ \pn\ 1 
Mat IOFat 
2 X 10-5 


Tr— m ^ 5 X 10^ 
\Bm2\ ) 


m ' 


10 - 


\pn\ 

Mn \B, 




m ' 


10-4 
\Bmjl"" 


m. 


(51a) 

(51b) 


IF I 

\Bf,Nj 


^ M]sr 
P and 


In estimate (51a I we assumed inequality (49), and in estimate (51b) we took into account relation (50), as well 
as identity (15) for Fat; at the end, we assum ed that the produced on-shell neutrinos Nj are semirelativistic, i.e., 
1 GeV). Using estimate (51b) and recalling that IB^at^P is the largest among the mixings IBoat^P, 

^hat for Mat = 1-5 GeV we can take IB^aTj p = |Bt-tv^. p, whose upper 
This implies that, at present, we can expect the values Lose 0.1-10 


|B,-aTj P, we can see from Table Fl 
bounds are given in the right column of Table^ 
m for the oscillation length. Of course, implicitly we assumed that the energies of the {B or B^) mesons, which decay, 
are not very high so that the assumption |pAr| Mjq would be justified. If |Losc| > 10 m, we would need quite a large 
detector, cf. Eqs. (47) and (48). 


4 Fat = (i/2)(rATi -t according to definition ( |l5[ l. 

^ We do not want to have \y\ > 10^ because the CF violation effects are then very suppressed m- 
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If we have \Lo 


0.1-1 m idet), our formulas (42)-(43) for LNV decays and (|45|)-(46l for LNC decays indicate 


that such oscillations can be detected and measured, once a sufficient number of such decays is detected, with the 
first (production) and the second (decay) vertices being within the detector. In this way, the oscillation length 
Lose oc I/AMat could be determined, and thus the mass difference AMjv (^ M^)- 

It is also interesting that these formulas indicate that in such a case the phases and could be measured 

as well. These phases could be determined, for example, by comparing the modulation of the measured differential 
effective decay widths \ L)/dL for the and B~ decays into /xctt, because the phase difference between 

the two oscillatory modulations is 2 x 02i\ cf. Eq. (42) for the LNV and Eq. (45) for the LNC case. The factor sin021 


appears in the CP asymmetry factor .Acp oc sin02i for these processes. For example, this asymmetry for the LNV 
case is 




^ e 7r+) — r(L+ —)■ /x+e+TT ) 


T{B- fj,-e-Tr+) + r(L+ /x+e+TT”) 


oc Psin0(f^)^ 


-f 1’ 


(52a) 

(52b) 


where y = AM^/Tn [cf. notation and factor P ~ 1 depends principally on the ratios of mixings \B(N 2 \/\BeNi\ 

{i = /X, e) and ratio /C1//C2 [cf. notation (18)]. This factor Acp can be substantial if y = AM^/^n is not too small 
in absolute value, e.g. if |y| ^ 10. We refer to Refs. [TSl EH] for more details on this. An interesting aspect here is 
that, by the described measurement of the angle 02i we could conclude that the CP asymmetry Acp is nonzero even 
in the case when jyj ^ 1, i.e., when this asymmetry is practically unmeasurable. 



FIG. 6: Graphical schematical representation of the differential decay rate —>■ yen-, L)/dL, cf. Eqs. (42) and (45l. 


The differential decay width dr^g^\B —>• /xe7r;L)/dL of Eqs. 
where L is the distance between the two vertices and Lmax = Ldet ■ 


(42) and (45) is presented schematically in Fig. p] 
n order to interpret how to measure this differentim 


decay width, we recall that this quantity is the limit (1/AL) x rgg'^^(P —yetr; L < L < L + AL) when AL —4-0 
(i.e., AL <C |Losc|), and here L is the distance between the production {y.-Nj) and the decay (Nj-e-n) vertex. To 
measure such a quantity, a sufficiently high number of events for each chosen bin L < L' < L + AL would have to be 
measured (with AL | Lose I and L < L^et)- 

There may exist another complication in such measurements. Namely the length Lose can vary in the detected 
events of the considered decays because Lose oc Pn'Jn oc \pn\ = \Pe+PTv\- In principle, the 3-momentum pN = pe+p-^r 
can be measured in each such decay, i.e., Lose can be determined in each such event. The graphical representation 
Fig. [prefers to a class of events which, among themselves, have approximately equal value of Lose, i-C-, approximately 
equal \pn\- If the decaying B^ (or B^) mesons were at rest in the lab frame, then the value of \pn\ is such a frame 
would be fixed by kinematics, namely 


\P_ 


4°) I = 


Ml 

MM 


Ml 

Ml 


(53) 
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where notation (6b) is used. 


In reality the B mesons coming into the detector have energies Eb > Mb- Let us assume that the incoming B 
mesons in the lab frame have all approximately the same 3-momentum pb = \'Pb\z parallel to the direction z of the 
tube of the detector where both vertices are detected, and that the detector tube is relatively narrow. Then the vector 
Pn in the detected events is the 3-momentum which can be obtained from = \pf^\z [cf. Eq. (53)] by a constant 
boost in the direction —z, bringing us from the B rest frame into the lab frame where B’s have the (approximately) 
constant 3-momentum |p_b|z. Thus the lab 3-momentum pjq = |p)v|i is approximately constant also in such a case. 
In such a case Lose would be approximately the same for all the detected events B —>■ pN —>■ peir in the tube, and the 
oscillation modulation indicated in Fig. could be measured, including the phase O 21 relevant to CP violation. 


V. DETERMINATION OF THE HEAVY-LIGHT MIXING COEFFICIENTS 

Measurement of the differential effective decay widths dT^°^^\B^ \ L)/dL can also lead to determination of the 
absolute values and of the heavy-light mixing coefficients. For example, for determination of I 

{j = 1,2) it is convenient to consider the decay widths for the semileptonic LNV decays B^ —>■ (and/or for 

the LNC variant) 




^ r{B+^ p+N)T{N ^ p+TT-) 


InPn 

y.[\B. 


+ \B^n,\^ + 2 \B^nA^\B^n,\^cos ± , (54) 


where = 2arg(i3^7V2)—2arg(i3^7Vi), in analogy with the/re case Eq. (14). Expression (54) differs from expression 

(42) only by the replacements p (and 621 ^ 4 ‘ 2 i)- In the case of B^-^ p^p^ir^, thesymmetry factor (1/2!) 
due to two identical muons in the final state gets canceled by the factor 2 coming from the inclusion of the square 


of the crossed channel amplitude, cf. Ref. [13]. The crossed channel amplitude in the decay B^ —>■ p^e^ir^ did not 
enter because it represents a different (distinguishable) process once the two vertices are identified and localized in the 
experiment. Stated otherwise, the processes B^ —)■ p^N —)■ p^e^ir^ and B^ —>■ —)■ e^p^ir^ are distinguishable 

once the vertices are identified. 

Measurement of the average (over L) /dL) and of the modulation amplitude A(dT^^'^'^/dL) determines the 


quantities (see also Fig. 


\B^NA^ + \B^N,f = {F), 2\B, 


= AF, 


(55) 


respectively. Let us denote as Ni the neutrino with larger mixing element (jS^ArJ > iF^Wal).® Then the heavy-light 
mixing coefficients (j = 1,2) are determined as well 


|2 _ 1 
fJ-Ni I — 2 

2 _ 1 
2 


\B, 

\BfiN 2 


+ AF+ yJ{F) - Af) , 
+ AF - a/(F) - Af) , 


(56a) 

(56b) 


VI. CONCLUSIONS 


In this work we considered the phenomenon of neutrino oscillations in semileptonic decays of B mesons via on-shell 
heavy nearly mass-degenerate Majorana neutrinos Nj {j = 1,2): the lepton number violating (LNV) decays B^ —>■ 
p^Nj —)■ p^e^TT^, and the lepton number conserving (LNC) decays B^ —)■ p^Nj —>■ Since the neutrinos 

contributing to such decays have to be on shell (the off-shell neutrinos give completely negligible contributions), the 
relevant flavor analogs are not and [Eq. (Q], but the truncated combinations Afi and A// [Eqs. ([^ and ([^], 
which are combinations of only the heavy mass neutrinos Ni and A2. The central results of the work are Eqs. (42) and 
(45) for the LNV and LNC differential effective decay rates dT^^‘^\L)/dL. These quantities must be interpreted as 


We have the freedom to do that, because the formulas are symmetric under the exchange A^i N 2 (i-e., (p 2 i —4>2i\ AM/v ^ —AMjy, 

-^osc ' ^ -^osc)* 
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(1/AL) X (L < L' < L + AL), where L' is the measured distance between the production vertex {p.-N) and the 
decay vertex (A-e-7r), and AL is considerably smaller than the oscillation length |Losc| = ‘2-tt\pn\/ Here, 
Pjv is the (approximately constant) 3-momentum of the intermediate A(,’s, and mass quantities are AMpf = 
where |AMw| ^ = M^- We argued that it is possible to have |Losc| 0.1-10 m if the 3-momenta ps and 

thus pn are not too large. If the detector length is comparable with |Losc|, and a sufficient number of mentioned 
decays is detected, we argued that it will be conceivable to measure the L-dependence of the differential decay width 
/dL, i.e., the oscillation modulation effects. By measuring these effects, the value of Lose could be discerned 
and thus the value of the mass difference AMpf. Moreover, by measuring such effects it would be possible to discern 
the phase 6*21 [cf. Eq. (14) and (19)], which plays an important role in the CP violation. In addition, magnitudes 
iL^jVil and \BppfJ (£ = p,, e) of the heavy-light mixing coefficients could be measured. 

In all the presented formulas, we can replace the initial meson by any other heavy pseudoscalar meson (such 
as B^, Df, B^), and final meson by any other lighter meson (such as , Df). This is performed by simply 

replacing everywhere the meson masses, the decay constants and the CKM elements accordingly [cf. Eqs.([^, etc.]. 
Among the initial mesons, those which can be copiously produced are evidently preferred. Additionally, those with 
higher mass are preferred because then the masses of the on-shell neutrinos can be larger; thus, the probability 
for the decay within the detector, Pn{L), can be more significant. Also, the preferred initial mesons are those which 
have less CKM suppression, i.e., whose CKM element |Vq^q^| in Eq. (6a) is not too small. Therefore, the preferred 
initial mesons are, in general, Df (IKcsj ~ 1 and Md = 1.97 GeVjand B^ i\Vcb\ ~ 0.04 and Mb = 6.28 GeV), 
but not necessarily ([Kt,] « 0.004 and Mb = 5.28 GeV). 

Eurthermore, as shown in Sec. |Vj the obtained formulas can be extended in a straightforward way to the decays in 
which the two charged leptons are equal, i.e., M^ —)■ —)■ and M^ 


p^N, 
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Appendix A: The quantum mechanics approach to oscillation 


In this appendix, we show that the amplitude approach to on-shell oscillations in the considered processes, as 
presented in the main text of this work and following mainly the amplitude approach of Ref. |36j , is consistent with 
the usual (quantum mechanics) approach to neutrino oscillation [37] (cf. also [3B]) applied to these processes (within 
the approximations used in such approaches). 

We recall that the relevant e- and /i-flavor analogs in the considered processes are the combinations ( |3^ of only 
the two almost mass-degenerate heavy neutrino eigenfields Nj (j = 1,2), because the other components (including 
the light neutrino mass eigenfields 1 ^ 2 , 1^3) are off shell or are assumed to be off shell in the considered processes. 
Following the usual (quantum mechanics) approaches to neutrino oscillation, cf. [37] (cf. also [3^ 1. the e and p heavy 
flavor analogs Afa (a = 1,2) of the heavy neutrino mass eigenstates Nj (j = 1,2) [cf. Eq. ([^] are represented as 
quantum mechanical states [cf. Eq. (3a) for the corresponding fields] 


K) = bI,\N,) + bI^\N2), 

Fa) = B^,\N,) + B^2\N2) (a = 1,2), 


(Ala) 

(Alb) 


where in Eq. (Alb) we assume that the physical neutrinos Nj are Majorana. Here, we use notation (3b) for the 2x2 


matrix B with normalized lines. In the wave function approach [37], these wave functions are in the Schrodinger 
representation and, consequently, have the following evolution in time t: 


Weit)) 

Fa(t)) 


2 2 2 


2 


2 2 


exp(-iL;jt)| A,) = exp{-iEj 

i=i /3=it=i 


t) IV,), 


(A2a) 

(A2b) 
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where we recall notation (3b) used for the 2x2 matrix B, and the inverse matrix is, consequently, 


B ^ 

DetJS 


B 22 —^12 

—-S 21 -Sii 


(A3) 


and B is the complex conjugate of this. In Eq. (A2) the notation Ej = E^. is used for the energy of the neutrino 
mass eigenstate \Nj), where E^^ is given in Eq. (29). I’he states |fVj) (j = 1,2) are orthogonal to each other 


{N,\Nk) = S,k. 


(A4) 


We note that the 2x2 matrix B matrix, Eq. (3bI, although having its two lines normalized, is, in general, not unitary, 
and therefore 


{JV, |A4) = ^ ^ 0, {JV, IAA 2 ) = ^ B,,B2, + 0, 


(A5) 


i.e., the states of th e he avy flavor analogs, |A/’a) and/or lAC/?), are, in general, not mutually orthogonal. As a 
consequence of Eqs. (A4) and these flavor analogs are normalized states 


(AAilAAi) = (AA2IAA2) = 1 = (AAilAAi) = (AA2IAA2). 


(A6) 


In the LNC decay —>■ /r+e“7r+. Fig. the neutrino flavor state produced in the first (production) vertex is jA//), 

and the state disappearing at the second^ecay) vertex is |A/i), cf. Fig.[^ Therefore, the relevant oscillation amplitude 
in this decay is {Afi\Af 2 {t)). 



FIG. 7: The LNC decay —>■ tt"*": at the production vertex, the |A/i) state is produced; at the decay vertex, the IA/ 2 ) 

state is absorbed. 


Using relations (A2a), (A5) and (A6), we obtain^ the following expression for the relevant oscillation amplitude 

(AAi|A4(t)): 


(AAi|A 4 (t)) = {eM-^EN^t)Bl,[(^B* 

. B11B21 + B12B22 


+ exp{-tEN,t)B;2 [{b* + (U* ( 

_^ ^ 

exp{-iEN-,t)B2i ^B22 - B^2 [ 


1 

DetS* 

+ exp{—iEj^^t)B 22 


B 11 B 21 


B 12 B 22 


-B 


21 


B 


11 


B 11 B 21 + B 12 B 22 


(A7a) 


(A7b) 


Transforming Eq. (A7al into Eq. (A7b), we use for B the complex conjugate of identity (A3). In this quantum 
mechanics approach, the terms in Eq. (|A7|) wi th ex p(—iEf^.t) correspond to the terms exp(—’ z) of the corre¬ 
sponding amplitude A{B~^ —>■ /r+e“7r+) m Eq. (44a). If the two approaches are to be consistent with each other, then 
the ratio of the coefficients at exp{—iEMA) exp{—iENA) in Efl- (A7b) is equal to the ratio of the coefficients at 


^ The algebra is performed in analogy with the usual quantum mecha nics approach to light neutrino oscillations [37], except that now 
we have, in general, the nonorthogonality of the two flavor states Eq. (A5 I. 
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exp(—ipATj • z) and exp(—ip]Y^ ■ z) in Eq. (44a|. This means that, for consistency, we need to have® 

0 


^22 ^12 


BuB 


11-021 


B12B22 


-Bh+Bl^ (S11B21+B12B2 


Bii 

Bi2 


(A8) 


By direct cross-multiplication, it is straightforward to check that this identi ty r eally holds. In checking this identity, 
it is enough to use only the normalization of the lines of the B matrix, Eq. (3b): |i?iiP -I- I.B 12 P = 1. 

In an analogous way, we can check that this quantum mechanics approach is consistent with the amplitude approach 
of the main text in the other cases: 

• In the LNC case B~ —» M~ e~*~7r~: in the explanation above (Eig. [^, the states IA/ 2 ) and |A/i) get replaced by 
|A/’ 2 ) and lA/"!); cf. Eq. (A1). 

• In the LNV case B+ —)■ /i+e+7r“: in the explanation above (Eig. [^, the state |A/i) gets replaced by lA/*!); 
cf. Eq. ( |Al| ). 

• In the LNV case B~ ^“e“7r+: in the explanation above (Fig. [^, the state IA/ 2 ) gets replaced by |A/’ 2 ); 
cf. Eq. (lATl). 
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